First-principles study of polarization in Zni-^Mg^O 
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Wurtzite ZnO can be substituted with up to ~30% MgO to form a metastable Zni-^Mg^O alloy 
while still retaining the wurtzite structure. Because this alloy has a larger band gap than pure ZnO, 
Zni-^Mg^O/ZnO quantum wells and superlattices are of interest as candidates for applications in 
optoelectronic and electronic devices. Here, we report the results of an ab-initio study of the sponta- 
neous polarization of Zni-^Mg^O alloys as a function of their composition. We perform calculations 
of the crystal structure based on density-functional theory in the local-density approximation, and 
the polarization is calculated using the Berry-phase approach. We decompose the changes in polar- 
ization into purely electronic, lattice-displacement mediated, and strain mediated components, and 
quantify the relative importance of these contributions. We consider both free-stress and epitaxial- 
strain elastic boundary conditions, and show that our results can be fairly well reproduced by a 
simple model in which the piezoelectric response of pure ZnO is used to estimate the polarization 
change of the Zni-^Mg^O alloy induced by epitaxial strain. 

PACS numbers: 77.22.Ej, 77.65. Bn, 77.84.Bw 



I. INTRODUCTION 

Recently, much attention has been paid to wurtzite 
Zni^Mg^O alloys as candidates for applications in op- 
toelectronic devices in the blue and ultraviolet region. 
ZnO is a wide-band-gap semiconductor with a direct gap 
of ^3.3 eV. The band gap becomes even larger if Zn 
atoms are substituted by Mg atoms, which have a similar 
ionic radius, allowing the construction of quantum- well 
and superlattice devicesi Similar behavior is well known 
for the zincblende GaAs/Al^Gai-^As system and is the 
basis of much of modern optoelectronics^ Recent trends 
have led in the direction of fabricating similar structures 
in wide-gap semiconductor systems such as wurtzite III- 
V nitrides^ and in Zni-^Mg^Oi^ There has also been 
recent interest in other kinds of nanostructures based on 
the ZnO and Zni-^Mg^O materials systems 

Pure ZnO prefers the wurtzite crystal structure, while 
MgO adopts the cubic rocksalt structure. Substitution of 
Zn by Mg results in a metastable wurtzite alloy for cer- 
tain magnesium concentrations. Experimental reports 
concerning the growth of these alloys on sapphire sub- 
strates indicate that Mg concentrations up to ~30%, 1,5 
or even ~50%jili can be achieved. 

Many ab-initio calculations of the properties of 
the parent compounds MgO and ZnO have appeared 
in the Iiterature>ii*i2iiii4 The properties of ternary 
Zni-^Mg^O alloys have been less well studied. There 
have been calculations of the dependence of the band 
structure and band gap on concentration xi^ Regard- 
ing the question of crystal structure and stability, Kim 
et al. has shown that the wurtzite Zni-^Mg^O alloy is 
stable with respect to the corresponding rocksalt alloy 
for x < 0.375>i& Similar results were obtained by Sanati 
et al. but for x < 0.33r 7 However, Sanati et al. also have 
shown that Zni-^Mg^O is unstable with respect to phase 
separation into wurtzite ZnO and rocksalt MgO phases 
even for low x values. This means that Zni-^Mg^O al- 



loys are not thermodynamically stable, consistent with a 
rather low observed solid solubility limit for Mg in ZnO>i£ 
The success in fabricating samples with higher concen- 
trations indicates that the phase separation is kinetically 
limited, i.e., the time scale required for the alloy to phase 
segregate into the two lower-energy constituents is long 
compared to the growth time at the growth temperature. 

To our knowledge, there have not been any previ- 
ous calculations of the the polarization properties in the 
Zni_a;Mg x O system. This is an important property to 
study, since if an interface occurs between a ZnO re- 
gion and a Zni^^Mg^O region within a superlattice or 
quantum-well structure, bound charges are expected to 
appear at the interface. These charges, in turn, will 
create electric fields that are likely to affect the electri- 
cal and optical properties of the quantum- well devices. 
In the present work, therefore, we have undertaken a 
study of the polarization and piezoelectric properties of 
Zni_ a; Mg x O. 

The structure of the paper is as follows. In the next 
section we describe the computational methods used in 
our work. In Scc. lIIII we introduce the six supercell struc- 
tures that were constructed and used as the structural 
models for the alloys of interest. Then, in Sec. IIVI we 
report the main results of this work. Finally, a brief 
summary is given in Sec. IVl 



II. COMPUTATIONAL METHODS 

Calculations of structural and polarization properties 
are carried out using a plane-wave pseudopotential ap- 
proach to density- functional theory (DFT). We use the 
ABINIT code package^ with the local-density approxi- 
mation (LDA) implemented using the Teter parametriza- 
tion of the exchange-correlations* and with Troullier- 
Martins pseudopotentialsiSi For the Zn pseudopotential 
the 3d valence electrons are included in the valence, as 
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their presence has a significant effect on the accuracy of 
results, 22 A plane-wave basis set with an energy cutoff of 
120 Ry is used to expand the electronic wave functions. 
A 6 x 6 x 4 Brillouin-zone /c-point sampling is used for 
pure wurtzite ZnO, and equivalent /c-point meshes are 
constructed for use in all wurtzite supercell calculations. 
The electric polarization is calculated using the Berry- 
phase approach^ 



III. SUPERCELL STRUCTURES 

In the present work we study the properties of six dif- 
ferent models of the ternary Zni-^Mg^O alloy, to be de- 
scribed shortly. However, first consider pure wurtzite 
ZnO. It can be viewed as two identical hexagonal closed- 
packed {hep) lattices; we take the O sublattice to be 
shifted in the +z direction relative to the Zn sublattice. 
Three parameters determine this structure: a and c are 
the lattice constants of the hep lattice, and u describes 
the shift between the two sublattices. 

Replacing some of the Zn atoms by Mg atoms, we get 
a ternary Zni-^Mg^O alloy. Of course, the real alloy 
is highly disordered. In order to carry out calculations 
using periodic boundary conditions, we construct ordered 
supercells having the same Mg concentration x as the 
alloy of interest. By comparing properties of different 
supercells having the same x, we may obtain a rough 
estimate of the size of the errors that result from the 
replacement of the true disordered alloy by an idealized 
supercell model. 

When constructing supercells, we restricted ourselves 
to structures having hexagonal symmetry about the z- 
axis, since real Zn^^Mg^O alloys have this symmetry 
on average. This makes the calculation and interpreta- 
tion of the results easier. We constructed six model alloy 
structures: one for x — 1/6 (Model 1), two for x = 1/4 
(Models 2 and 3), one for x = 1/3 (Model 4) and two for 
x = 1/2 (Models 5 and 6), as follows. 

The simplest alloy one can make (Model 5) is obtained 
by replacing the Zn atoms by Mg atoms in every second 
Zn layer along z, giving a structure with Mg concentra- 
tion x = 1/2 and retaining the primitive periodicity of 
pure ZnO (four atoms per cell). Similarly, if one replaces 
every fourth layer of Zn by Mg, one arrives a model with 
x = 1/4 (Model 2); this has an eight-atom supercell with 
the primitive lxl in-plane periodicity but with a doubled 
periodicity along the z-direction. 

In the remaining models, we retain the primitive peri- 
odicity along z but expand the size of the supercell in the 
x-y plane, as illustrated in Fig. ^ Models having 2x2 
in-plane periodicity (Models 3 and 6), and those having 
\/3 x V3 periodicity (Models 1 and 4), are specified with 
reference to Figs. Qfa) and (b), respectively. Models 3 
and 6 thus have 16 atoms per supercell, while Models 1 
and 4 have 12 atoms. In Model 3 we assign b = c = Zn 
and a = d = Mg in Fig. Qfa), obtaining a model with 
x = 1/4 in each cation layer and x = 1/4 overall. Model 




FIG. 1: (Color online.) Top view of cation layers of super- 
cell models for Zni-^Mg^O alloys, (a) Structures with 2x2 
periodicity (Models 3 and 6). (b) Structures with v3 x Vo 
periodicity (Models 1 and 4). Atoms 'a' and 'b' lie in the top 
cation layer, while 'c' and 'd' are one layer below (see text). 



6 corresponds to b = d = Zn and a — c = Mg; this results 
in alternating cation layers with x = 1/4 and x = 3/4, 
for an overall Mg concentration of x = 1/2. Turning to 
the V3 x v3 structures in Fig. IHb), one can see that 
the hexagonal symmetry requires that all atoms must be 
the same (c atoms) in one of the layers. We construct 
Model 1 by assigning b = c = Zn and a = Mg, yielding 
alternating layers with x = and x = 1/3 for an average 
x = 1/6. Finally, for Model 4 we set a = c = Zn and 
b = Mg so that the layer concentrations are x = and 
x = 2/3, averaging to x = 1/3. 

Of course, it would be possible to generate more su- 
percell models of the alloy by expanding the periodicity 
or reducing the symmetry. However, the six models de- 
scribed above provide a reasonable coverage of concen- 
trations in the range < x < 1/2 with some redundancy 
(for x = 1/4 and x — 1/2). We have thus chosen to limit 
ourselves to these six models in the present work. 



IV. RESULTS 

A. Pure ZnO and MgO 

To determine the crystal structures and cell parame- 
ters of pure ZnO and MgO, we carried out DFT calcula- 
tions for both materials in both the wurtzite and rocksalt 
structures. For wurtzite ZnO we obtained lattice param- 
eters a = 3.199 A, c = 5.167 A and u = 0.379. While 
these results are very close to previously reported the- 
oretical values^ they slightly differ from experimental 
values^ (a = 3.258 A, c = 5.220 A and u = 0.382). The 
cohesive energy (defined as the energy per formula unit 
needed to separate the crystal into atoms) is found to be 
8.26 eV. Comparing this to the cohesive energy of rock- 
salt ZnO (8.03 eV), one may conclude that ZnO prefers 
the wurtzite structure, in agreement with experiment. 
For rocksalt MgO we found a = 4.240 A and a cohe- 
sive energy of 10.00 eV. We find that if we start with a 
plausible wurtzite MgO structure with a, c and u simi- 
lar to that of ZnO, the crystal can monotonically lower 
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TABLE I: Theoretical equilibrium lattice parameters for bulk 
ZnO and for models of Zni-^Mg^O. Subscript 'free' indicates 
zero-stress elastic boundary conditions, while 'epit' indicates 
that a is constrained to be identical to that of bulk ZnO (the 
values in column V are thus identical by construction). 





X 


«free (A) 


( C / a )free 


a cp it (A) 


(<=/a) ep , 


ZnO 


0.0 


3.199 


1.615 


3.199 


1.615 


Model 1 


0.17 


3.216 


1.605 


3.199 


1.624 


Model 2 


0.25 


3.230 


1.593 


3.199 


1.625 


Model 3 


0.25 


3.225 


1.600 


3.199 


1.628 


Model 4 


0.33 


3.238 


1.589 


3.199 


1.630 


Model 5 


0.5 


3.266 


1.564 


3.199 


1.635 


Model 6 


0.5 


3.256 


1.580 


3.199 


1.640 



TABLE II: Theoretical cohesive and formation energies (eV 
per formula unit) for bulk ZnO and MgO and for each super- 
cell model. 



ZnO O0 8^258 O0 - 

Model 1 0.17 8.496 -0.053 

Model 2 0.25 8.602 -0.093 

Model 3 0.25 8.612 -0.083 

Model 4 0.33 8.729 -0.123 

Model 5 0.5 8.955 -0.176 

Model 6 0.5 8.958 -0.173 

MgO 1.0 10.004 0.0 



its energy along a transformation path in which a in- 
creases, c decreases, and u tends toward 1/2 in agree- 
ment with the previous results of Ref. 0. The minimum 
occurs at u = 1/2, which corresponds to the higher- 
symmetry h-MgO structured For this structure we ob- 
tain a = 3.527 A and c = 4.213 A, in good agreemenliiiii^ 
with previous calculations. We find its cohesive energy 
to be 9.81 eV, consistent with the fact that MgO prefers 
the rocksalt structure. (For more details concerning the 
previous theoretical literature on lattice parameters and 
binding energies, see Ref. HH l 

The main goal of the present work is to study the po- 
larization and piezoelectric properties of Zni-^Mg^O. 
For reference, our calculated spontaneous polarization 
for pure ZnO is found to be — 0.0322 C/m 2 , and its 
piezoelectric coefficients are e 31 = —0.634 C/m 2 and 
e 33 = 1.271 C/m 2 . Note that the value of the sponta- 
neous polarization differs somewhat from the previous 
theory of Dal Corso et al.yO* who reported a polarization 
of —0.05 C/m 2 when using the experimental u — 0.382; 
our value becomes much closer to theirs if we also use the 
experimental u. Since we are primarily interested in dif- 
ferences of the polarization with respect to pure ZnO, we 
do not believe that these small discrepancies are impor- 
tant. The values of piezoelectric coefficients are in good 
agreement with previous theoretical calculations of Wu et 
alM who found e 3 i = -0.67 C/m 2 and e 33 = 1.28 C/m 2 
(and who also provide comparisons with other theoretical 
and experimental results). 



B. Crystal structure and energies of alloys 

For each model described in Sec. IIIII we calculated the 
hep lattice parameters a and c in the equilibrium state. 
Since we are interested in properties of Zni-^Mg^O lay- 
ers that might be grown on a ZnO substrate, we also 
calculated the lattice parameters for epitaxially strained 
structures (i.e., a fixed to that of pure ZnO). The results 
are given m Table IJ In both cases, the c/a ratio exhibits 
an almost linear dependence on x. However, this ratio 
is found to decrease with x for the fully relaxed struc- 



tures, while it increases with x when the epitaxial strain 
condition is enforced. 

In Table ^ we give cohesive and formation energies 
for each alloy. One can see that in every case the for- 
mation energy is negative. Thus, according to our LDA 
calculations, at zero temperature the Zni-^Mg^O alloy 
is never stable with respect to phase-separated wurtzite 
ZnO and rocksalt MgO. (Of course, at T > a small 
solid solubility of Mg in wurtzite ZnO is expected^) 

C. Polarization and piezoelectric properties 

The results of the calculations of spontaneous polariza- 
tion are given in Tabic ITTT1 both for the fully relaxed and 
for the epitaxially strained cases. Note that the values 
of polarization for models having the same x are fairly 
consistent with one another; the choice of supercell does 
not significantly affect the overall trend with x, which is 
reasonably smooth. A linear fit P(x) — P(ZnO) + Ax 
yields coefficients of A frcc = (-0.088 ± 0.009) C/m 2 and 
A opit = (0.024 ± 0.002) C/m 2 . The latter value may be 
of direct interest for experimental studies of epitaxial su- 
perlattices and quantum wells. 

Thus, with increasing Mg concentration x, the abso- 
lute value of the polarization increases for the relaxed 
structures and decreases for the epitaxial structures with 
fixed a. This behavior is very similar to what we saw 
in Sec. HVfjl for the c/a ratios, suggesting that the c/a 
ratio may be a dominant factor in determining the total 
polarization. Indeed, since 2e 3 i + e 33 ~ 0, one expects 
the polarization to be almost independent of a change 
in volume (isotropic strain), so that the change of c/a 
should be the most important strain effect. 

In order to study more thoroughly the role of strain 
and other factors in determining the polarizations of the 
Zni^Mg^O structures, we first define AP t ot to be the 
polarization of the alloy superlattice structure relative to 
that of pure ZnO. We then decompose APtot into "elec- 
tronic," "ionic," and "piezoelectric" contributions as fol- 
lows. First, we construct an artificial Zni-^Mg^O super- 
lattice structure in which the structural paramters (a, 
c, and all internal coordinates) are frozen to be those 
of pure ZnO, and define AP e i ec to be the polarization 
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TABLE III: Calculated values of total polarizations of 
Zni^Mg^O alloy models (C/m 2 ). Subscript 'free' indicates 
zero-stress elastic boundary conditions, while 'epit' indicates 
that a is constrained to be identical to that of bulk ZnO. 
Superscript 'est' indicates value estimated by the model of 
Eq. 1). 





X 


-ffrcc 


Pcpit 


pest 
* cpit 


ZnO 


0.0 


-0.0322 


-0.0322 




Model 1 


0.17 


-0.0423 


-0.0277 


-0.0279 


Model 2 


0.25 


-0.0501 


-0.0247 


-0.0247 


Model 3 


0.25 


-0.0470 


-0.0244 


-0.0250 


Model 4 


0.33 


-0.0565 


-0.0230 


-0.0239 


Model 5 


0.5 


-0.0789 


-0.0199 


-0.0222 


Model 6 


0.5 


-0.0699 


-0.0202 


-0.0225 



of this structure relative to that of pure ZnO. Next, we 
allow only the internal coordinates of the Zni-^Mg^O su- 
percell to relax, while continuing to keep a and c frozen 
at the pure-ZnO values, and let APi on be the polariza- 
tion change produced by this internal relaxation. Finally, 
we allow the lattice constants to relax as well, and de- 
fine APpiczo to be the associated change in polarization. 
Clearly AP = AP olcc + AP ion + AP piczo . 

The results of such a decomposition are given in Table 
II VI for the stress-free case. For scale, recall that these are 
changes relative to P(ZnO)=-0.0322 C/m 2 . The purely 
electronic contributions AP c i 0C are quite small, showing 
a relatively poor correlation with x. The contribution 
APi on associated with the ionic relaxations is also quite 
small, although it is typically 2-3 times larger than AP i oc 
and shows a clearer trend (becoming more negative with 
increasing x) . By far the largest contribution comes from 
the piezoelectric effect of the strain relaxation, being typ- 
ically 5-10 times larger than the ionic one. A similar table 
can be constructed for the case of epitaxial strain; its hrst 
four columns would be identical to Table llVl because of 
the way AP on and AP c i cc are defined, and the values in 
the remaining columns can be deduced from the infor- 
mation given in Tables |ml and |IV| The results indicate 
that the piezoelectric contribution also dominates in the 
epitaxial-strain case. 

This being the case, it seems likely that many of the 
polarization-related properties of the Zni_ x Mg x O alloy 
can be estimated by using a model based on the piezo- 
electric effect alone. For example, one might hope that 
SP = P C pit — Pfroo, the difference between the epitaxially- 
constrained and free-stress polarizations at a given x, 
could be estimated by a linear approximation of the form 

r ,-> rt a cpit ^free , Ccpit Cf ree , * 

6P = 2e 3 i he 33 — . (1) 

In fact, we find that this is the case even if we use the 
piezoelectric constants of bulk ZnO, already obtained in 
Sec. IIV Al in this formula. Using the computed value of 
Pfi-cc reported in the third column of Table IIIII together 
with the constrained a values and epitaxially-relaxed c 



TABLE IV: Theoretical values of electronic, ionic, piezoelec- 
tric and total contributions to polarization (C/m 2 ) for each 
model, relative to bulk ZnO. 





X 




APion 


A-fpiczo 


APtot 


ZnO 


0.0 


0.0 


0.0 


0.0 


0.0 


Model 1 


0.17 


0.0001 


-0.0022 


-0.0081 


-0.0101 


Model 2 


0.25 


0.0018 


-0.0023 


-0.0175 


-0.0180 


Model 3 


0.25 


0.0000 


-0.0027 


-0.0122 


-0.0148 


Model 4 


0.33 


0.0009 


-0.0038 


-0.0214 


-0.0243 


Model 5 


0.5 


0.0023 


-0.0063 


-0.0427 


-0.0467 


Model 6 


0.5 


-0.0019 


-0.0062 


-0.0296 


-0.0377 



values given in the last two columns of Tabled we report 
the computed estimates P^pit = Pf r0 o + SP in the last 
column of Table IIIII The use of the piezoelectric coef- 
ficients of pure ZnO is not obviously justified except at 
small x, but the results show excellent agreement with 
the computed P ep it values in the fourth column even up 
to x = 0.5, where the error is only about 10%. This 
approximation thus seems to work quite well. 



V. SUMMARY 

We have investigated the polarization-related prop- 
erties of wurtzite Zni^Mg^O alloys using calculations 
based on density-functional theory in the local-density 
approximation and the Berry-phase approach to calcu- 
lating electric polarization. In particular, we have stud- 
ied the dependence of the spontaneous polarization on 
Mg concentration using six alloy supercell models with 
hexagonal symmetry, spanning the range of Mg concen- 
tration from x — 1/6 to 1/2. We performed these cal- 
culations both for free-stress and epitaxial-strain elastic 
boundary conditions. 

Our results indicate a roughly linear dependence 
of spontaneous polarization on Mg concentration, al- 
though the sign of the linear coefficient is opposite 
in the free-stress and epitaxial-strain cases. In or- 
der to understand this behavior in more detail, we de- 
composed the change in polarization into electronic, 
lattice-displacemcnt-mcdiatcd, and strain-mediated com- 
ponents, and found that the latter component is dom- 
inant. This means that the change in polarization is 
mostly governed by piezoelectric effects connected with 
the x-dependent changes of the a and c lattice constants. 
We further confirmed this picture by showing that the 
polarization changes could be well approximated by a 
model in which the only first-principles inputs to the 
model are the piezoelectric coefficients of pure ZnO and 
the ^-dependence of the equilibrium lattice constants of 
the Zni-^Mg^O alloy. These results suggest that charg- 
ing effects associated with polarization discontinuities 
in ZnO/Zni-^Mg^O superlattices and quantum wells 
should be subject to prediction and interpretation in a 
fairly straightforward manner. 
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